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ABST'EAC^ 


Data for the system n-pentane/u-heptane on porous Chrom sorb-102 
adsorbent -were obtained at 150, 175, and 200° C for mixtures containing 
zero to 100£ n-pentane by weight. Prior results showing limitations on 
superposition of pure component data to predict ml ti cccponcat chroma- 


tograms vere verified. The tbemodycgsic parameter nR^ 


Kas found to 


be a linear function of samp le composition. A noaporous adsorbent 
failed to separate the system because of large input samp le dispersions. 
A proposed automated data processing scheme involving magnetic tape 
recording of the detector signals and processing by a mini computer was 
rejected because of resolution limitations of the available A'D conver- 
ters. Preliminary data an porosity and pore size distributions of the 
adsorbents vere obtained. 
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Tho work done for this report is part of a long term effort 
at Rensselaer (l, 2 , 3 , 4, 5 , 6) to develop fundamental design 
techniques for a gas chromatograph . Specifically, the overall 
project is to develop methods for designing a. -id optimizing a 
gas chromatograph-mass spectrometer system for use in chemical 
analysis to be carried out in an unmanned roving vehicle on the 
surface of Ears, although the results are by no means limited to 
this application. 

A gas chromatograph is a chemical separation device which 
utilizes the phenomena of adsorption-desorption. Chemical 
species, carried by an inert or carrier gas, successively adsorb 
on and desorb, from a packed bed of granular particles (with or 
without a liquid substrate) at different rates. Because of this, 
a multicomponent sample injected into a chromatograph will elute 
as a series of separate chemical species. 

In the present steadies, a packed column is operated isothormally 
with a constant flew rate of carrier gas (helium). A sample is 
injected into the carrier gas and as it passes into the chromato- 
graph, the various species diffuse, adsorb, and desorb. For each 

• • • * » 

• 

chemical, axial diffusion is represented by the Peclet number, 
mass transfer to the adsorbent surface is represented by the 
parameter N^qq, and adsorption-desorption is represented by the ' 
thermodynamic parameter xeRq . The Peclet number is a function of 
both the fluid mechanics and physical properties of the system; 
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®tOQ primarily a function of the fluid mechanics; and mHo I s a 
• function of the chemical species and the -column used* 
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An existing gas chromatograph test facility has been operated 
iso thermally to generate data on both single component and binary . 
mixtures of n-pentane and n-heptane. A previously developed 
chromatographic data reduction computer program and an equilibrium 
adsorption model were used to analyze the data. 

Binary data obtained on a porous column shows that superposition 
of- predicted single component behavior is only a first order 
representation of actual binary data. The parameter ieRq is found 
to bo a linear function of composition. This suggests that the 
adsorption follows a non-linear isotherm, and that there are 
possible interference effects in binary systems. 

Data obtained on a non-porous column is well predicted by 
the model, and fails to show spreading that was observed on 

porous columns. This shows that the mechanism of intraparticle ! 

* 

diffusion is responsible f^r the observed spreading. 

- The porous packing, 0-102, was found to consist of porous 
microspheres. The manufacturer reports an averc -e micropore 
diameter of JO A, while experimental results show the macropores 
range from 400-700 A. 

* i 

A system of automatic data processing was evaluated and j 

i 

found to offer less accuracy than possible with the present \ 

i 

manual processing. Minor changes in the detector control j 

circuit are recommended to improve manual data handling. - j 









PART III 
BACKDROBND 

Previous work on the gas chromatograph system Includes 
development of several mathematical models (1,2, 4, 6 ). Each 
• model was developed from the fundamental second-order partial 
differential equation set describing the chromatograph , and each 
predicts the behavior of a single component sample at the column 
outlet based on an impulse injection. Two second-order models 
which include z ial diffusion have been developed, the equilibrium 
adsorption model (2) and the non-equilibrium adsorption model (4,6). 
The former assumes pointwise adsorption equilibrium throughout 
the column (N^qq = CO ) , while the latter assumes a non- 
equilibrium situation (N^OCr = finite). Woodrcw (6) has shown that, 
for the conditions of the test facility used, the non-equilibrium 
model offers only minor improvement is prediction of output 
chromatograms. Because of this, and because the non-equilibrium 
model requires more extensive calculations, the equilibrium adsorp- 
tion model is used for data analysis in this report. This model 
assumes pointwise eq uili brium between the gas phase and the 
adsorbed phase of the packing, i, e. , the column is very long. 

. The equilibrium adsorption model is represented by 


-Pe(6-^) 2 /4f3 

which gives the response at the end of the column to a unit inpulse. 
To obtain experimental data, a previously built test facility ' 




was used ( 5 ). i‘ne data are analyzed using a general data 
reduction computer program (GIRP) prepared by Benoit (5)» which 
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given chromatographic input data, output data, and system 
parameters such as the Peclet number and will reduce the 

data to dimensionless form, predict an output chromatogram based 
on the model, and numerically and graphically compare the 
prediction with the actual data. This program predicts chromatograms 
for finite pulse injection by numerically convoluting the input 
data with the calculated impulse response. The physical parameters 
(Pe, Re, v, etc.} are Calculated by another computer 

program ( 6 ). 

Ir order to extend the model to multi-component systems, 
suj 3 imposition of single component chromatograms has been investigated 
by Keba ( 6 ). His work has shown that superposition offers only 
a first order approximation of actual binary data, and that 
parameter mR© is apparently composition dependent. Fur the rk 

has been done on binary systems in order to verify this result 
and to determine the variations of mRo with composition. 

The data of Benoit (5) and I'^ba ( 6 ) has shown spreading of the 
out-put chromatograms that is not predicted by the model. Since • 
a porous packing was used in obtaining this data, it has been 
postulated that this spreading is due to intraparticle diffusion(S) . 
In order to study this, data have been obtained on a non-porous 
column. 

Work is being done by Woodrow ( 7 ) to improve the mathematical 
model to obtain more accurate prediction. This requires data on 
the physical characteristics of the adsorbents, such as particle 
pore size distribution and column void fraction. Preliminary 





estimates of these properties were obtained from experimental 
work and from manufacturer ' s literature, 

> In an attempt to improve consistency in the data and to 
reduce the manual effort in handling the data, a method of 
automatic data process.^ was investigated, 

* x 
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RESULTS AND DXSCUSSICS 

A. SOPElPOSmCU AND HDJAHI DATA 

The mathematical model assuages a linear Isotherm for each 

component of the form (l) 

. ■ 

• 

y - . 

This Is generally a valid assumption if the sample composition in 
the carrier gas is small. For a typical sample used in this work, 

J) 

the sample compositions are on the order of 10~ mole fraction 
(see Appendix A). 

This isotherm implies that equilibrium for each component is 
independent of any other components present. Because of this, it 
seems reasonable to use a linear superposition of the single 
component predictions to predict the behavior of a binary mixture. 

Data on the binary system of n-pentane and n- heptane were ob- 
tained on the porous Chromosorb 102 column (C-102) described in 
Table I, at 150°C, 175°W and 2Gj°C. The sample compositions 
ranged from an excess of pentane to an excess of heptane. In 
all cases a syringe was used for injection of two microliter 
liquid samples. The experimental conditions for the binary data 
are given in Tablo U. 

Analysis of binary data shows that the tins of appearance of 
the output peaks in binary mixtures lags the time predicted from, 
superposition of pure component predictions. This effect is 

i 

shown for a typical experiment an Figure 1. This time lag 



TABLE I 


COLUMN CHARACTERISTICS ’ 


- 

Chromosorb IQ? 

Carbouax 1500 

Length (cn) 

100 

100 

Outside Diamster{cn) 

0.32 

0.32 

Inside Diarcoter(cn) 

0.22 

0.22 

Hash 

. 60 /S 0 

fo/80 

particle Size(cm) 

0 . 025 / 0.0177 

0 . 025 / 0.0177 

Compos i tion 

>H_croporous 
Styrene-Divinyl 
Benzene Polymers 

Polyethylene Glycol 
20 p by weight 

Tempera tore Range 

Room-250°C 

Roosn-lSOOC 

Application 

Separation of lew 
molecular weight, 
very polar sub- 
stances such as 
water , alcohols , 

Separation of high 
boiling polar com- 
pounds. General 
purpose column for 
polar samples. 


amines, acids, 
esters, ketones, 
ethers, aldehydes, 
and to a limited 
extent, light gases. 



TAELS II 


EXPrJUKBJTAL (XKDITKKS 
fcr bhuki d\ta cu c-ia2 


Temperature 
deg C 

Column Inlet 
Pressure, psig 

Helium Flowrate 
cu cm/rdn 

- 150 

30.0 

■ : 29.4 

175 

30.0 

30.0 

200 

30.0 

26.4 


* at 70°P and lU.7 psia 


occured for both components for all mixtures and at all temperatures 
It was found that the time lag increased as the component 
composition became more dilute. (see Table ill) 

The peak tin® is primarily related to the thermodynamic 
parameter eRq (2). This parameter cannot be estimated a priori 
and is determined from actual data by a curve fitting technique 
which varies e & 0 until the peak time of the predicted chroma to— • 
gram matches the peak time of the output data (5). To determine 
eRq values for binary chromatograms, the output of each component 
is analyzed separately by use of the GUP. By use of this method, 
mR Q values have been calculated for each component. As a general 
trend, the value, of mRo decreases as the component becomes more 
dilute at constant temperature. From Figures 2 and 3. eRq is 
seen to be a U mar function of composition. It should be noted 
that some binary samples ( 755 “ and 995 “ pentane by weight) wore not 
suf' > ' uiently separated at 200°C on C-102, so thRq values could not 
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TABLE HI 

• TDE LAGS EJ HIXAHI DATA 

Pentane-Eeptane System on C-1Q2 
Flow Rate of 2g.4 cc/min at 150 °C 


Comoosltion 


Time of Pentane 


Time of Heptane 


weight Jo 

Peak, sec 

Peak, sec 

100 H 


285.5 

-99 H - IP 

95.0 

294.0 

75 H - 25 P 

88.5 

307.5 

50 H - 9 ) P 

E3.0 

327-5 

25 H - 75 P 

84.0 

30+.O 

1 H - 99 P * 

81.5 

432.5 

100 P 

8 O .5 



. 



jfc. _J tL . 



Fig -ure 2 Effect of Composite 
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be calculated for these samples , 

Since Keba (6) also noted the occurrence of a time lag, eRq 
values were estimated for his data on C-102 at 150°C and 200°C 
by means of the following relation (5) 

* 1 

. ZutvQ (X ' ■ 

®max,o ~ ®max,i “ * 

t 

which gives a value within ± 10 ^ of the actual value* These values 
appear in Figures U and 5, and it can be seen that the same general 
trends are exhibited. 

There are several possible reasons for the observed behavior 
of the parameter mRo. First, the observed variations with 
composition are probably due to changes in the equilibrium constant 
m,' since Rq is essentially constant at constant temperature and 
pressure. This indicates either that the isotherm is not ? -near 
as assumed but there is 1 no interf erence effect, or that the 
isotherm is non-linear and dependent upon the concentrations of 
both components. 

The first explanation makes use of the fact that, since the . 
binary samples were all two microliters, the- composition of a 
component in the carrier gas decreases as its composition in the 
sample decreases, A typical non-linear non-interactive isotherm 


is 


*L,1‘ 


Ky? 

1 + Kyj; 


a Langmuir, single-component isotherm. This isotherm would give 


flEtETC* 



Figure 4 Effect of Composition on Heptane mR, 
• for Pentane-Heptane System on C-102 
(Data of Keba) 




Figure 5 Effect of Compos ition on Pentane iriR, 
for Pontane-Heptano System on C-102 
(Data of Keba) 
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a different value of m for each value of 7*, but does not include 
any effect of another component. If the adsorption followed this 
isotherm, changes in the pure component injection volume would 
produce variations in the value of m. Specifically, if y* is 
reduced, m will decrease, which is consistent with the observed 
effect. If this is indeed the case, the above equilibrium relation 
can be used in the mathematical model, and superposition can be 
used to predict binary data, although its use is not rigorously 
correct. Use of this non-linear isotherm would considerably 
complicate the calculations. It may not be possible to obtain an 
analytical solution, and numerical methods may be required. 

Any interference effect in the binary systems is probably in 
the fora of competition between the two components for adsorption 
sites. This can be represented -by a non-linear interactive isotherm 



C L,1 


%yf 

1 + K iy j . + K 2 y|' 



a binary Langmuir isotherm. The calculations for this non-linear 
interactive model are not much more complicated than those for the 
non-linear, non-interactive model. Use of this model, if it is 
correct, would predict binary chromatograms and the assumption of 
superposition would not be required. 

B. TEMPERATURE EFFECTS 

In general, the value of mR 0 increases as the temperature 
increases. The form of this temperature variation has been 
proposed as (6) 







nR 0 oC exp(-E/RT) . l 

Correlation of pure component C-.102 data according to this ! 

expression is given in Figure 6, Over the limited temperature 
range considered, the data does follow the correlation and 

activation energies on C-102 are on the order of 11 kcal/gm-mole j 

[ 

for pentane and 15 'kcal/gm-mole for heptane. These values indicate . I 

■ ’ i 

i 

that the controlling mechanism is probably the physical adsorption i 

process. Pentane is well modeled even at 150 °C (Figure 7)* but 

the quality of prediction for the heptane data increases at 

higher temperatures as shown in Figures 8, 9» and 10 (see Appendix 

B for computer output) . This substantiates the fact that physical 

adsorption is controlling. ’ 


C. PRESSURE DEPENDENCE OF mR 0 • > 

The parameter mR 0 consists of two separate factors: m, an 

equilibrium adsorption constant; and R 0 , the ratio of moles of "j 

gas in the bed to moles of liquid in the bed. ' . • > 

The parameter m appears as a result of the assumption of a \ 

•. s 

i 

linear isotherm. This is expressed as i 

i 

■ x L = Kp* , ^ j 

but assuming an ideal gas (8) one obtains ■ ' 

xl - KPy* 
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Effect ‘of Temperature on 
Parameter mR 0 for C-102 
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Figure 9 Prodicted and Actual 'Chromato gratis 
■ for Haptano on C-102 at I 75 °C 


























If X Is assumed to be Independent of pressure, u Is seen to be 
inversely proportional to pressure, , 

The parameter Rq is defined as (l) 

R o = €fr/F(l -€)g L - 

Far an ideal gas, the density is proportional to pressure (at 
constant temperature) . The void fraction of the bed is independent, 
of pressure. F is the stationary phase volume which is liquid. 

Both this term and the. liquid density are very weak functions of 
pressure. 

Combining these effects, one obtains 

nfio- me^/F^I-^^^Ci/pKp) <p° . 

Therefore, as a rough estimate, one would expect mHo to be 
independent of pressure. This is borne out by the data obtained 
(see Figure 6).. . - 

D. DATA CN HCH-PCROUS PACKING - 

A Car tow ax l^OO column (C-l^OO) described in Table I, which 
is believed to approximate a n on-porous packing., was used to 
generate single component data for pentane, heptane, and acetone. 
The experimental conditions for this data are given in Table IV. 

It is seen in Figures 11 and 12 that the model very closely 
predicts the shape of the actual data. This indicates that the 
presumed, mechanism of intraparticle diffusion is probably 
responsible for the observed spreading of the porous column data. 



V 





TAELS IV 

EXPERIMENTAL CamiTICKS 
FCR DATA CSJ C-I5OO 

Temperature Column Inlet Helium Flowrate* 

deg C ' Pressure, psig cu cm/min 

.100 20.0 31,6 

125 • 20.0 . 29.6 

150 20.0 27.7 

* at 70°^ and 14.7 psia 

Analysis of these data presented Sousa problems. Pentane and 
heptane are weakly adsorbed on C-I5OO, resulting in a very sharp 
impulse response. However, there is much tailing in both the 
input and output data so a relatively coarse time step was used in 
the computer analysis. The combination of these two factors 
accounts for the failure of the convolved response to predict 
exactly the height of the actual data in the tao figures (11 and 
12 ) shown. The actual area under the convolved response is 
somewhat less than that under the actual chromatogram. Since the 
area under the curve is proportional to the amount of the sample 
injected, both of these areas should be the same for proper 

i 

prediction. Because of this, each value of the convolved 
response should be multiplied by a constant so that the area 
under this curve is identical- to the area under the actual data. 

It can be seen that if this were done, the two curves would be 
virtually identical. 




dieted a 1 
■ Heptane 















■ . j. — - A 

28 . 

V 

The observed oscillations at the end of the convolved response 
(Figure 12) are probably due to the large time step used and 
errors in the data (see following section). These oscillations 

are inconsequential in the determination of mRo and the predicted 

; 

response. 

Over the range of temperature studied, the C-I 5 OO column 
would not separate the binary mixtures of pentane-heptane. In 
order to obtain a qualitative estimation of the degree of 
separation obtained, moment analysis was used. Voytus (2) has 
shown that the time of the center of gravity of a chromatogram (mean) , 
as well as the variance, can be obtained by moment analysis from 
the Laplace transform solution of the model and knowledge of the 
system transport parameters. The mean is defined as 

/'ey d0 

Q = . . 

r ° yde 

- •'o 

> . 

The variance, or second moment about the mean, is 

0-a = L (e-efy £ e - • 

: c *. ie 

For the equilibrium adsorption model, it lias boon shown (2) 

that 


and 


0 = ^ = 1 + lMo 

= p a (2/Pe) 


f 






whore is the variance of the impulse response. The variance 
. of an actual chromatogram (9) is 

where q£ is the variance of the input pulse. The value of 
is obtained from acuual input data using the defining equation 

A 

given above. Values of C£ for the input data of Keba have been 
calculated ( 7 ) and are used in this analysis. Voytus (2) has 
shpwn that at least 94 ^ of the total area of the chromatogram is 
within a 20 spread about the mean. 

The analysis, done here is for the pentane-heptane system at 
100°C, Pure component values of eRq were used and the binary 
system was represented by superposition of the two pure component 
predictions. Table V summarises the values of the parameters used. 


TABLE V 

PARAMETER VALUES ECR 
MOMENT ANALYSIS 



rrR 0 

Pe 

a ?, 

oi. 

CJ- , 

3 *1 

cr , 

cr, sac. 

Pentane 

0.4566 

IO7&O 

19.63 

0.03 

19.66 

8.88 

Heptane 

0.3315 

12500 

19.63 

o.o4 

19.67 

8.88 


The chromatograms were assumed to be approximately symmetrical 
about the mean ( implying that the mode is about equal to the mean). 
With these assumptions, the semi-qualitative Figure 1J was drawn 
(see Appendix C for calculations). It shows the groat degree of 
overlap of the pentane and heptane chromatograms , and that, while 
some separation occurs, it would not be detectable. 



Prediction 








V 


, i v 




31 


The. temperature behavior of the C-1500 data is shown in 
Figure- l4, which presents only a quantitative representation. It 
is obvious that the straight lines drawn are subject to a large 
degree of error, since there are so few data points. The v>lue 
of mR 0 for acetone and heptane increases as temperature increases, 
while the' value for. pentane remains relatively constant over the 
temperature range studied. The kRq values for pentane and heptane 
diverge as the temperature decreses, indicating that a dectec table . 
separation might occur at a temperature below 100°C. The 
activation energies for pentane, heptane, and acetone on both 
C-1500 and C-102 are given in Table VI (see Appendix D for 
calculations). The energies are much lower for the C-1500 column 
than for the C-102. This is because the components are very weakly 
adsorbed on the C-I5OO. This could also explain the erratic behavior 
of the pentane and heptane mRo. ‘ \ 

_ TABLE -VI . 

• ACTIVATICW ENERGIES 

■ ■ kcal/gm-mole . . 


C-1500 

Pentane 0. 

Heptane * 


C-102 


11 . 


6 ) 


Acetone 


0.5 

2.0 ' 


15 

6 


(Ref. 



Figure 1 H Effect of Temperature on Parameter mR 0 for C-I50O 



£. SOURCES OF EERQR. 

Since the data are presented in a form which requires 
computer calculations, there are many sources of error. Errors 
are associated with, reading the value of the chromatogram from the 
strip chart (see Part IV, 0), assuming isothermal column operation 
determining helium flew rate, estimating the various parameters 
(Pe, Re, etc.), using numerical methods in computer calcula- 

tions, and with computer roundoff. When the parameter mRc, is 
calculated, all these errors' are present. Since mR 0 is determined 
by a curve fitting technique, the major source of error ^ in the 
processing of the actual data. 

Each data point has two dimensions associated with it: 
composition and time. The height of the trace on the strip chart 
above the base line is proportional to the composition (y) of the. 
component in the carrier gas, • This height is measured using an 
engineer's rule so the error here is +0.02 inch. Since the 
maximum deflection obtained was about 3 inches for most of the 
experiments , the average error is about 1$, For the data studied 
in thus report, a value of y is read every half second. Since 
time lines appear only once a second on the strip chart, there is 
an error in estimating the time of about +0.1 second. Another 
error associated with time is the location of ’'zero 0 time. This 
is generally taken to be the time of appearance of the input 
pulse, but since values are taken only once every half second, 
there can be an error of as much as 0.4 second. The manner in 
which theso errors are propagated in the calculations is difficult 


to determine, although it seems likely that the value of eRq 
obtained is only accurate to within + 20p. 

Possible ways to reduce this error are in obtaining a full 
scale deflection on the strip chart (see Part IV, 0) and using 
a time line of 0.1 second. The latter is possible if the paper 
speed is increased to the next faster speed. 1 in/sec. This 
reduces the error, in determining the time associated with a 
point, but uses the light sensitive paper at five times the current 
rate.' As this paper is expensive, the improvement brought about 
by such a change should be evaluated carefully to determine if it 
is justified. 

F. PACmO CHARACTERISTICS 

An effort was made to obtain no , , iical characteristics of 
the C-102 and C-I500 packings used, Inforiiatiun on composition 
and particle size appear in Table I. 

It was found that C-102 is manufactured by the Rohm and Haas 
Company under the trade name u Amberlite XAD-2" and that Perkin- _ 
Elmer merely dries and crushes it before packing the columns (10). 

A technical b ull etin and samp", e of XAD-2 were obtained from Rohm 
and Haas, and a sample of C-102 was obtained from Per-' cin-Elner . 

•The structure of XAD-2 appears in Figure 15, while soma of its 
properties are listed in Table VII, From the concept of Figure 
15, it is postulated that each particle consists of a large 
number of very small microspheres, which, in turn, are porous. 

The C-I5OO packing consists of firebrick (manufactured by 
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TABLE VH 

HiOPERTIES 0? XAD-2 
(Reference 11 ) 


Appearance 

Hard, Spherical 
opaque beads. 

• 

Solids 

51 to 55 

PorositjKml pore /ml bead - dry basis) 

0.40 to 0.45 

Surface Area(m^/go - dry basis) 

330 

Effective Size (mm) 

. 0.30 to 0.45 

Harmonic Dean Particle Size (mm) 

. 0.45 to 0.60 

Average Pore Diameter (A - dry basis) 

90 ’ 

True Wet Density in distilled 
water (gra/ml) 

1.02 

Skeletal Density (gm/ml) 

i.o; 

Bulk Density (lbs/ft^) 
(gm/cc) 

4 o to 44 
0.64 to 0.70 

Bulk Density*, (gm/ml) 

0.4o 


* experimentally determined. It is assumed that the bulk density 
reported by Rohm and Haas is the density as shipped, i.e. , about 
• 5 o£ wet. 




Union Carbide) doped with polyethylene glycol by Perkln-ELmer (10). 
A sample of this packing was also obtained. 

A mercury penetration porosineter was used to estimate 
actual pore distribution and external void fraction of these samples 
It is probable, however, that the results obtained are subject to 
large, errors, as the proper sample holder for these powder samples 
was not available. The available sample holder, designed for . 
large pieces of material, was used and the results are considered 
to be semi-quantitative. The results of the tests for XAD-2 
and for .C-102 appear in Figures l 6 and 17 respectively. Meaningful 
results were not obtained with the C -1500 sample, because some of 
the sample flowed out of the holder during the experiment. 

In Figure 16 , the section of the curve up xo about 100 psia 
represents the external void volume, while the region between 
1000 psia and 10,000 psia represents the macropores within the 
particles (most probably the volume between the mi exospheres) . 

From these results, the mean diameter of these macropores is 
approximately 0.0^ micron (4O0A). Further, since the penetration 
Volume is still increasing with pressure, it appears that more 
pores still exist within the sample ( 10,000 psia is the limit of 
the apparatus). 

Inspection of Figure 17 shows that the C-102 sample exhibits 
the same characteristics as the XAH-2, indicating that the two 
materials are probably the same. The region between 1000 psia 
and 4000 psia represents the macropores, which have a mean diameter 
of about 0.07 micron (700 A), in fairly good agreement with the 
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XAD-2 data. The sharp change In slope between 8000 psia and 
10,000 psia indicates the presence or more pores of scalier 
dimension. 

Basically, the only difference in the two curves is in the 
region below 100 psia. This section represents the external void . . 
volume for both samples. In the XAD-2, this volume is practically 
filled even at the lowest pressure ( 3 psia) , while in the C-102 
greater pressure is required to fill the voids. This is because the 
particles of C-102 are much smaller than the XAD-2 (about l/20 
the size) , so that the space between C-102 particles is much 
smaller. 

Rohm and Haas reports a mean pore diameter of 90 A. It is 
possible that this is the mean diameter of the micropores (i.e., 
the pores of the microspheres). Assuming a contact angle of l^O 0 , 
a pressure of about 20,000 psia would be needed to fill these pores 
in the porosimetcr. Although this pressure cannot be reached with 
the apparatus, the presence of these micropores is indicated by 
the XAD-2 and C-102 data. • . . 

Based on the pores ime ter results, the void fraction and 
specific surface area of C-102 were calculated (see Appendix E) 
to be approximately 0.33 and m 2 /gm respectively. The value of 
surface area is much less than the reported value of XOO-UOO m 2 /gm 
(10). This demonstrates two things: 1) the existence of pores 
which would be filled at pressures above 10,000 psia; and 2) the 
large error associated with these results because of the sample 
holder used. Because of this, the results are considered as a 
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seni-quantitative representation of the actual properties. 

The parosiraeter is a. valuable tool in determining the 
properties of the packings. With the proper saaple holder, 
accurate measurements of pore distribution, external void fraction, • 
surface area, and density can be obtained. 

G. AUTOMATED DATA PROCESSED . • ' ' ' ' 

In an attempt to improve consistency in data analysis and to 
reduce the manual effort in data handling, the automatic data 
processing scheme shown in Figure IS was evaluated. This figure . 
also shews the currently used “manual" method of data processing. 

In the manual data processing method, the electrical data 
from the composition detectors at both the column entrance and 
exit are recorded on a strip chart by means of a vis ic order 
oscillograph. This machine also prints a time scale. Data 
points are then manually selected from the strip chart at 
specific time intervals, using an engineer's rule. For data 
taken in this work, a 0.5 second time interval was used. The . _ 
discrete data must then be punched on computer cards for final 
analysis using the computer. . 

In the proposed automated processing method, the electrical ■ 
data would be recorded both on magnetic tape by an FK recorder 
and on the strip chart. This is done for visual certification of . 
the data. The FK recorder would then be taken to a minicomputer 
and the data read into the machine through an analog-to-digital 
converter. This produces discrete data points with the specified 
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time interval between points. The minicomputer would then 
produce both a printed listing and a punched tape of .the data 
points. The punched tape would be read directly into the main 
computer for final data analysis. 

The detector control circuit appears in Figure 19, showing 
the addition of connectors used for the automated system. Since 
neither side of the Wheatstone bridge is at ground potential, . -d 
the FJI recorder available (see Appendix F.for specifications) 
recces voltages with respect to ground, two differential 
amplifiers were required. As none were available,, an analog 
computer was used in this capacity for purposes of system evalu- 
ation. To develop specifications for the differential amplifiers, 
measurements of the voltage across the two sides of the bridge 
were made under various conditions. The results of these tests 
appear in Table VIII. In order to obtain full output ( 5 volts) 
from the Fl-I recorder, differential amplifiers with gains of 1000 
are required for samples with dilute components, while gains of 
100 will suffice for ether samples. 

The A/D converter of the minicomputer available for use has 
a range of -10 v to +10 v and reads voltages at increments of 
78 mv„ This allows a signal to be broken down into 256 parts at 
most. Since the output of the FH recorder is only 0 v to 5 v, 
any data recorded could be resolved into only 64 parts. Resolution 
of 256 parts could be obtained by biasing the 5X recorder output 
with 2.5 vdo and using a four gain amplifier to provide a signal 
which ranges from -10 v to +10v. 





CARLE DETECTOR CONTROL CIRCUIT 





V 



-ii: 


table vxrr 

i 

VOLTAGE OUTPUT FROM 
WHEATSTONE BRIDGE 


‘C-102 COLUMN 


Pulse . 

Voltage, mv 

Input* (maximum) 

30.- 45. ; 

■ Output (maximum) 


. Pure Heptane 

15. 

Dilute Heptane 

0.45 

50 $ Heptane 

10. 

Pure Pentane 

30. 

Dilute Pentane 

0.35 

50 $ Pentane • 

15. 

c-isoo column 


■ _ Pulse 

Voltage, mv 

Input* (maximum) 

45. 

Output* ( maximum) 

4o. 

NOISE LEVELS 


Output Channel 

+ 0.02 mv 

Input Channel 

+ 0.01 m 

* approximately the same for all 

compononvs 




A .chromatogram on an oscillograph strip chart with full 
scale deflection ( 5' inches) .can be read with an engineer's 

rule to the nearest one-hundredth of an inch, giving a. 
maximum resolution of about 500 parts. This is about twice the 
maximum resolution obtainable from the automatic processing 
system. It should be noted that the automatic processing will 
provide consistent .data reduction from year to year. However, 
this consistency can be obtained at the price of : l) about $1000- 
$1500 for two required ‘differential amplifiers; and 2) less 
accurate data reduction. It is therefore recommended that the 
proposed automatic data processing s.etem not be implemented and 
that attention be directed to obtaining full scale deflection on 
the oscillograph strip charts. 

This can be accomplished by providing variable attenuation 
in 1 ■ control circuit (Figure 19). For most samples, an 
attenuation of one causes the trace to go off scale, while an 
attenuation of two gives less than full, scale deflection. Use 
of a variable potentiometer in conduction with the present ten-' • 
position potentiometer should provide the capability of • 
obtaining full scale deflection for almost every sample. 

Since the manual processing system is time consuming and 
subject to errors (see Part IV , E) , other systems for automatic 
data processing should bo explored. Use of a digital voltmeter 
with printing capabilities seems feasible, although the co3t is 
high. Such a system should receive more careful attention in 


future work 



CONCLUSION AND RSCOJENDATI 


1. Linear superposition offers only a first order represen- 
tation of actual binary data. The parameter mR 0 is a linear 
function of component composition. 

2. It appears the actual isotherm is non-linear. Studies 
should be made of the effect of sample size on mR 0 to determine 
if the isotherm is non-interactive or interactive. Once this is 
determined, the correct isotherm should be used in the mot el and 
changes made to give binary predictions. 

3. For pentane and heptane, mR 0 is proportional to exp(-E/RT) 
in the temperature range considered. On the porous column, 

these data are better :'odeled at higher temperatures. 

4. The parameter rcR 0 is pressure independent. 

5. Heptane data on a non-porous column are well modeled, 
showing that intraparticle diffusion is probably responsible for 
observed spreading of heptane data on r. porous column. 

6 . The C-I 5 OO column would not separate the pentane-heptane 
system at the temperatures considered. Operation of this column 
at lower temperatures may give separation. Another non-porous 
column should bo considered. 

7 . The C-102 and XAD-2 appear to be the same material. 

Each particle is composed of porous macrosphores. The tnacropores 
are on the order of 400-700 A, while the micropores have a mean 
diameter of 90 A. A sample holder for powders should be obtained 
and the porosimeter should be used to determine the packing 
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character'll ics . 

8. The proposed automatic data processing system should not 
be implemented, as the cost involved is not justified. 

9. Bettor resolution of the data can be obtained by obtaining 
full seal© deflection of the trace on the strip chart. This can 
be done by .using variable potentiometers in the control circuit. 

10. Many errors are associate^ v?ith data processing, 
especially in manual handling of the data. The value of mR 0 is 
assumed to be accurate within 20;1. 

11. Other systems of automatic data processing should be 
investigated, and the manual system should be evaluated to 
eliminate sources of error. 
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NCfEiCLATURE 


E- 

¥ 

L 

m 

bRq 


H 


tOGr 


Pe' 


P 

P 

R 

Re 

t . 
T ■ 


v 


X L ' 

y 

y* 


activation energy, kcal/gm-mole 

fraction of stationary phase volume which is liquid 
length of column, cm 

adsorption equilibrium constant = y*/* L 
thermodyimraic parameter ~ * 

number of transfer units, dimensionless 
Peclet number, dimensionless ' . 

partial pressure . 

total pressure 

gas constant, 1.987 .cal/gm-mole 
Reynolds number, dimensionless 

ratio of moles of gas to moles of adsorbate in bed. 
time, sec 

temperature, deg K 
carrier gas velocity, co/sec 
concentration in adsorbent phase 
concentration in vapor phase 
concentration in vapor phase at equilibrium 


GREEK LITERS 

\ 

p = i +i/na 0 

€L - bed porosity 

^ - density 

O' - - standard deviation 

- dimensionless time = vt/L 
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^aax,o ■ ’ - .dimensionless time of output peak 
e max,i - dimensionless time. of input peak 

SUBSCRIPTS 

0 - gas phase 

L — liquid phase 

J - unit impulse response 
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PART IX 

, * . l 

APPENDICES' 

A. Calculation of Sample Composition in Carrier Gas 
1 ■ ’ 

B. . Computer Output 

C. Calculations for Moment Analysis 

. D. Activation Energy Calculations 

E. Calculations for C-102 Surface Area and Void Fraction 
. F. Specifications .for FH Recorder 


The appendices appear in the reference; Meisch, A.J., 

"Gas Chromatography: Evaluation of Binary Data and 

Characterization of Adsorbent Properties," M. Eng. Project 
Report, Rensselaer Polytechnic Institute, Troy, New York, 
May 1973* Copies of the material may he obtained from 
Dr. P. K. Lasbmet, Systems Engineering Division, . 
Rensselaer Polytechnic Institute, Troy, New York, 12l8l. 
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